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Q . Abstract 

^ . There is a huge gap between properties of red-sequence selected massive galaxy clusters at 2: < 1 and 

^ I Lyman-break selected proto-clusters at z > 3. It is important to understand when and how the z > 3 

^— I. proto-clusters evolve into passive clusters at z < 1. 

We aim to fill this cluster desert by using the space-based NA^Ajim) imaging with the AKARI. The 

\ z' — N4 color is a powerful separator of cluster galaxies at z > 1, taking advantage of the 4000 A break and 

■ the 1.6/im bump. We carefully selected 16 promising cluster candidates at 0.9 < z < 1.7, which all show 
^vq I obvious over-density of galaxies and a prominent red-sequence. 

l/^ ■ At this redshift range, the mid-infrared Si^^rn/ Sg^rn fiux ratio is an extinction-free indicator of galaxy 

' star formation activity due to the redshifted PAH emission lines (6.2,7.7 and 8.6/im). We show statistically 

^5 . that the cluster galaxies have a lower Si^^rn/ S^^rn flux ratio than field galaxies, i.e., cluster galaxies already 

' have lower star-formation activity at 0.9 < z < 1.7, pushing the formation epoch of these galaxy clusters 

00 . to a higher redshift. 

■ Key words: galaxies: evolution, galaxies:formation, galaxiesiclusters: 
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\ 1. Introduction scale structure forms. High redshift clusters can help con- 
^ ■ strain feedback processes caused by star-formation and 
High redshift clusters are highly important objects in AGN (Silk & Rees 1998). 
variety of science: For example, it has been known that However, there has been a missing link at 1 < z < 2 in 
galaxy cluster environment strongly affects the evolution the current cluster surveys. Traditionally, in the nearby 
of galaxies. The morphology-density relation suggests the Universe, galaxy clusters are known to form a strong color- 
cluster environments transform galaxies into early- type magnitude relation (red sequence), dominated by red, 
(Goto et al. 2003b). The star formation rate of cluster early- type galaxies. And thus, galaxy clusters are pri- 
galaxies is significantly lower than in the field (Tanaka et marily found by searching for red sequence galaxies clus- 
al. 2004). A galaxy cluster is the most massive virialized tering on the sky (Gladders & Yee 2000; Goto et al. 2002). 
system in the Universe; cluster studies provide an impor- This method utilizes the 4000A break, and thus, optical 
tant tool for gauging the growth of structure and prob- surveys can only search for clusters at z<l (the limit set 
ing the density of the underlying dark matter and energy hy i — z color). This is why only few galaxy clusters are 
(Springel et al. 2005). A presence of (even a few) massive known at z > 1 today, most of them stemming from X-ray 
clusters at high redshift can constrain cosmological pa- surveys (Mullis et al. 2005; Stanford et al. 2005; Stanford 
rameters VLm and (Bahcall et al. 1997). Cluster studies et al. 2006). 

can illuminate how dark matter haloes collapse and large At z>3, Lya emission redshifts into the optical band. 

People found clustering of Lyman-a emitters (LAE) or 

Lyman Break Galaxies (LEG), often around a strong radio 
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galaxy or AGN (thus, these are biased surveys). However, 
these proto-clusters are groups of star-forming galaxies 
(LAE,LBG) and thus, fundamentally different from the 
low redshift clusters dominated by red early-type galaxies 
(Miley et al. 2004). It is important to investigate how 
and when these proto-clusters evolve into passive galaxy 
clusters in the nearby Universe. 

Also, due to these selection methods, there inevitably 
exists strong selection bias between the proto- and classic 
(passive) clusters. It is important to study the evolution 
of galaxy clusters using a sample selected in a consistent 
way. This paper provides one of the attempts to extend 
the classical cluster search to higher redshift. 

In summary, due to the red/blue wavelength limit of 
the optical surveys, there has been a hole in galaxy clus- 
ter surveys at 1< z <3. We aim to fill this missing link 
in previous galaxy cluster surveys by utilizing the deep 
near infrared imaging with the AKARI satellite com- 
bined with the ground based deep z'-band image with 
the SuprimeCam on board the Subaru telescope. Deep 
z'-hciiid images combined with the AKARI 4/im images 
straddle the 4000 A break of galaxies at 1 < z <2. And 
therefore, we can naturally extend the classical 4000A- 
break search of galaxy clusters to this redshift range of 
the cluster desert. 

Unless otherwise stated, we adopt the WMAP cosmol- 
ogy: {h.nrn^^L) = (0.7,0.3,0.7) (Komatsu et al. 2008). 

2. Data 

The AKARI, the Japanese infrared satellite (Murakami 
et al. 2007), performed a deep imaging in the North 
Ecliptic Region (NEP) in 2-24/im (14 pointings in each 
field over 0.4 deg^; Matsuhara et al. 2006; Matsuhara 
et al. 2007; Wada et al. 2008). Due to the solar syn- 
chronous orbit of the AKARI, the NEP is the only field 
with deep imaging in these wavelengths (but see Teplitz et 
al. 2005). Since the Spitzer space telescope does not have 
sensitivity at 9 < A < 20 jam (a gap between the IRAC 
and MIPS) except for the IRS peak up array at 16 jam 
and 22 //m, the AKARI NEP survey will be the only 
field with deep images in these wavelengths at least for 
a next decade. The depth provided by the AKARI satel- 
lite in A^4-band (4/im) reaches 8.0 /iJy in 5 (j(Wada et 
al. 2008), which can in principle detect bright galaxies 
up to z ^1.7. The 5 a sensitivity in other wavelengths 
{N2,N3,N4,S7,S9W,S11,L15,L18W and L24) are 14.2, 
11.0, 8.0, 48, 58, 71, 117, 121 and 275/iJy (Wada et al. 
2008). 

A part of the NEP deep field was observed in z'-band 
with the Subaru telescope (Imai et al. 2007; Wada et 
al. 2008), reaching the limiting magnitude of zab =26 
mag in one field of view of the Suprime-Cam. The 
combination of deep z- and A^4(4/im)-images provides 
us with a unique opportunity to search for galaxy clus- 
ters in the previous cluster desert at 0.9 < z < 1.7. The 
A^4(4/im)-band covers the same area observed with the 
N3,N4, S7, S9W, S'l 1 , L15, L18l^ and L24 bands, allowing 
us to investigate the mid-infrared properties of discovered 



clusters immediately. We use the AB magnitude system 
in this paper. 

3. Method 

3.1. Cluster Detection 

We have previously developed an efficient galaxy clus- 
ter finding algorithm called the Cut & Enhance method 
(Goto et al. 2002), which uses color-cuts and spatial clus- 
tering information to enhance a signal from galaxy clus- 
ters. In this paper we apply this method to z^ — colors 
of galaxies in the NEP deep field. 

In the right panels of Fig.l, we show the z' vs. z' — N4 
color-magnitude diagram. Orange dots are k-corrected, 
evolutionary models of a passive galaxies with IO^^Mq 
(total of the gas and stellar mass) as a function of redshift 
(indicated in the figure). This model assumes passive evo- 
lution of galaxies. 5 a detection limit of the N4 band is 
8.0 /iJy or 21.67 ABmag. This detection limit of the N4 
band is the deepest among the AKARI N2,N3 and N4 
bands. In addition, because of the redshifted 1.6/im bump 
(Sawicki 2002), the band does not become fainter as 
quickly as 7V2. As can be easily seen, colors of galaxies 
become redder and redder in z' — 7V4 color from z=0.9 
to z=1.7. As a reference, the overplotted dotted straight 
lines show the artificially redshifted color-magnitude rela- 
tion observed in the Coma cluster at z=0.9, 1.0, 1.1, 1.25, 
1.5 and 1.7. Because at z > 0.9, z'- and 7V4-bands strad- 
dle the 4000A break, the z^ — NA color becomes redder 
and redder with increasing redshift. In addition, the N4 
magnitude can take advantage of the 1.6/im bump and 
does not became as fainter as N2 or N3 bands at our tar- 
geted redshift range (Eisenhardt et al. 2008); in Fig. 2, we 
show three spectra of instantaneous burst model galaxies 
redshifted to z=1.3, with ages of 0.5, 1, and 10 gigayears 
of age with Salpeter initial mass function and the solar 
metallicity (Bruzual & Chariot 2003). The overlaid are 
the filter response of z\ N2,N3 and A^4 bands. At z=1.3, 
the redshifted 1.6/im bump is in the A/'4-band. This is why 
in Fig.l, the model track becomes redder, but not very 
fainter with increasing redshift (and thus moves straight 
up in the figure). Taking advantage of this color-shift, 
we can efficiently find galaxy clusters by using color-cut 
corresponding to each redshift, and by searching an over- 
density of galaxies within a color-cut (i.e., overdensity of 
galaxies with the same color). 

To find galaxy clusters in this field, we separate galaxies 
in the z' — N4 color space in the bins of A{z' — 7V4)=0.2. 
The width of this color-cut is set to take account of rel- 
atively large errors of the A^4 magnitude. We move this 
color-cut from blue {z' — A/'4=1.4 targeting z = 0.9) to red 
{z' — N 4=3.7 targeting 2; = 1.7) by 0.1 mag at a time as if 
scanning along the redshift. The adjacent bins have over- 
lap of 0.1 mag in color, assuring that a cluster which hap- 
pens to be on a border of a color cut is not missed. This 
color-cut removes most of fore- and background galaxies 
nicely, enhancing the density contrast of galaxy clusters. 

In each color-cut, we compute a local galaxy density 
of individual galaxies by computing a number of galax- 
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ies within 1.5 Mpc of radius within the same color-cut. If 
the galaxy has more than 2cr overdensity than the average 
galaxy density in the color-cut, the galaxy is selected to be 
in a dense region. We then connect these galaxies using a 
friend-of-friend algorithm with the linking length of (pro- 
jected) 0.5 Mpc (Goto 2005). If more than 4 galaxies are 
connected, we select them as a galaxy cluster candidate. 
Sometimes, a same cluster can be found in two adjacent 
color-cuts. In such cases, we combine such clusters if they 
were within 1.5 Mpc from each other and in the adjacent 
color-cuts. 

Using this method, we have found 16 galaxy cluster 
candidates in the range of 0.9 < z < 1.7. When tuned for 
lower redshift range, a known cluster, RXJ175719. 4+6633 
(z=0.6909; Gioia et al. 2003; Branchesi et al. 2006; Henry 
et al. 2006) was successfully detected with our cluster 
finder, showing our methodology works well. The list of 
these cluster candidates are presented in Table 1. 

3.2. Redshift Estimate 

The next important step is to measure redshift of each 
cluster we found. Spectroscopic redshift measurement is 
ideal. However, at 0.9 < z < 1.7, taking spectra of clus- 
ter galaxies which usually do not show emission lines is 
costly. Especially, for clusters at z > 1.5, near-infrared 
spectroscopy is required to measure Ca H&K absorption 
lines. For galaxies close to A/'4=22mag, such near-infrared 
spectroscopy is impossible with a reasonable amount of 
existing 8m class telescope time. 

However, cluster galaxies are known to have a tight 
color-magnitude relation (Kodama et al. 1998), i.e., 
galaxy clusters have many passive galaxies, forming a line 
in the color-magnitude diagram (called the red-sequence). 
Such red-sequence galaxies have a strong 4000A break, 
and thus, their photometric redshift can be measured 
much more reliably. Since galaxy clusters have many red- 
sequence galaxies, it is easier to measure photometric red- 
shift of galaxy cluster than those of individual galaxies. 
In the past, such cluster photometric redshift measure- 
ment reported accuracy of Sz < 0.02 using only 2 colors 
that straddle the 4000Abreak at the redshift (Gal et al. 
2000; Goto et al. 2002; Im et al. 2002). 

We apply this method to our clusters using the z' — 
color. First, we fit a line with the slope of the Coma color- 
magnitude relation to cluster member galaxies within the 
z' — N4: color-cut selected in Section 3.1. This model as- 
sumes the passive evolution of galaxies. Then, we re-select 
member galaxies within ±0.1 mag in z' — 7V4 along the fit- 
ted line, and within 1.5 Mpc from the cluster center, which 
are median coordinates of galaxies in the initial selection. 
We re-fit the line again to the new set of member galaxies 
to determine the photmetric redshift of each cluster. In 
most cases, a difference between initial redshift estimate 
and our final photometric redshift is less than 1%. 



4. Results 

4.1. RXJl 7571 9.4 + 6631 31 

In the AKARI NEP field, there is a known galaxy clus- 
ter, RXJ175719.4+663131, at z=0.6909. The ROSAT X- 
ray satellite has a similar sun synchronous orbit as the 
AKARI, and thus has deepest data around the NEP re- 
gion. This cluster was first found by the ROSAT satellite 
in X-ray (Gioia et al. 2003; Henry et al. 2006; Branchesi 
et al. 2006). The redshift of this cluster based on 5 spec- 
troscopic member galaxies is z=0.6909. 

In the upper panel of Fig. 3, we show z'- (left) and 7V4- 
(right) band images of this cluster. Galaxies marked with 
red circles have z' — NA color consistent to be at redshift of 
0.6909. It is seen that these galaxies are clustered within 
the 1.5 Mpc around the cluster center. In addition to 
galaxies marked with red circles, there are obvious over- 
density of galaxies in the z'-band image, showing this is 
truly a galaxy cluster. 

In Fig. 4, we show z' — NA vs z color-magnitude diagram 
of cluster galaxies with correct z' — NA color within the 
0.5 Mpc with large dots. A tight color-magnitude relation 
can be seen. The gray dots show all the galaxies within 
0.5 Mpc from the cluster center (Gioia et al. 2003). This 
cluster has also been successfully detected with our blind 
search for galaxy clusters. 

4-2. Most promising cluster candidates 

In this section, we present detected cluster candidates 
at 0.9 < z < 1.7. 

We have detected 22 cluster candidates in this redshift 
range in total using the method described in Section 3. 
Among them, four clusters are identical to the detection 
in the adjacent color-cut, and thus are merged into single 
clusters. Among the remaining 18 clusters, one was af- 
fected by an artificial spike of the A^4 image. Another at 
z=0.97 conincides with the position of the z=1.36 cluster, 
and can be considered as bluer galaxies of the background 
cluster. Therefore, the two candidates are removed from 
the list. In the left panels of Fig. 1, we show Subaru 
z'-band images of all 16 cluster candidates. The image 
scale is indicated in the green bar in the figure. Galaxies 
in red circles have di z' — NA color consistent with the red 
sequence at each redshift. Around these galaxies, there 
are obvious overdensity of galaxies including galaxies not 
detected with the AKARI, and thus, these candidates are 
likely to be real galaxy clusters. Some of cluster candi- 
dates even have a possible cD galaxy. In the central panels 
of Fig. 1, we show corresponding AKARI A^4 images of 
clusters. We show a corresponding color-magnitude rela- 
tion of these cluster candidates in the right panels of Fig.l. 
Filled circles are galaxies in S{z^ — NA) < 0.2 from the best- 
fit color of the red-sequence, and within 1.5 Mpc from the 
cluster center. The gray dots are all galaxies within 1.5 
Mpc around the cluster in the NEP deep SuprimeCam 
field. The orange dots are k-corrected, evolutionary mod- 
els of a passive galaxies with M = 10^^ (total of the 
gas and stellar mass) as a function of redshift (indicated 
in the figure). The model shows that NA magnitude does 
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Fig. 2. Spectra of galaxies with ages of 0.5, 1 and 10 gi- 
gayears redshifted to z=1.3. The spectra are from the in- 
stantaneous burst model of Bruzual &; Chariot (2003) using 
the solar metallicity and the Salpeter initial mass function. 
Overlaid are the filter response curve of z' ,N2,N3, and A^4. 

not become much fainter in this redshift range because of 
the 1.6/im bump. 

To quantify the significance of the overdensity of cluster 
candidates in Table 1, we have measured a local overden- 
sity, = cluster Infield ^ of clustcr corcs (galaxics within 
0.5 Mpc of radius around the median position of cluster 
galaxies). Here T^duster is a number density of cluster 
galaxies within 0.5 Mpc of radius and in the color-cut. 
S field is computed for each galaxies outside of the cluster 
in the same way as cluster galaxies (around 0.5 Mpc and in 
the color-cut), then we take mean to compute H/ieZd- The 
mean field galaxy density was computed using all galaxies 
except those belong to the cluster in the color-cut corre- 
sponding to the cluster redshift. The results are presented 
in Table 1. The mean 8y: is 5.5, which means the galaxy 
density of cluster cores is 5.5 times larger than that of the 
corresponding field galaxies. 

We also measure the significance of overdensity, a^] = 
[T^ciuster - field) /cf{T^ field)' We mcasurc fiuctuation of 
field in field)' The results are presented in Table 1. 
The mean cty: is 2.4, which statistically occurs with less 
than 2% of probability. We do not claim our sample is 
completely free from the contamination. True identity of 
the clusters needs to be verified spectroscopically. 

4.3. S15/S9 flux ratio 

In this subsection, we investigate star-formation activ- 
ity of cluster galaxies found in Section 4.2 using the PAH 
emissions redshifted into 15/im-band. The 5 a sensitivity 
in the AKARI S^W and L15 bands are 58 and 117 /iJy 
(Wada et al. 2008). With these sensitivities, the lower 
limit of the luminosity range we explore at 0.9 < z < 1.7 
are approximately > 2 - 9 x lO^^L© and > 1 — 7 x lO^^L© 



Fig. 4. The z' vs z' — color-magnitude relation of the 
RXJ175719. 4+663131. Galaxies plotted with the large dots 
have consistent z' — N4 color to be at z=0.6909 and within 
0.5 Mpc from the cluster center. The small dots are all galax- 
ies within within 0.5 Mpc from the cluster center. The red 
line is an expected location of the color-magnitude relation 
at z=0.6909. The orange dots show model galaxy colors with 
evolution as a function of redshift. 

in 9 and 15/im, respectively. 

In this redshift range of 0.9 < z < 1.7, star-forming galax- 
ies are known to have a large Si^^rn/ S^^rn flux ratio since 
prominent PAH emissions at 6.2,7.7, and S.Gjiim are red- 
shifted into the L15-band and form a pseudo-break to the 
continuum at < 5/im (see Takagi et al. 2007; Takagi et 
al. 2008 for a model prediction, and Ohyama et al. 2008 
in prep, for spectroscopic confirmation). And of course, 
the mid-infrared Si^^rn/ S^^rn flux ratio is free from dust 
extinction unlike optical. As a star formation indicator, 
we prefer the Si^^rn/ S^^rn flux ratio over the mere lu- 
minosity in L15, since 15/im luminosity largely depends 
on galaxy mass, and thus, is not an ideal tracer of the 
star formation activity when targetting massive galaxies 
in clusters (Goto et al. 2002; Goto 2005). In the pres- 
ence of warm dust of the AGN, it is difficult to tell if the 
L15 luminosity is from the AGN or the star- formation. 
On the other hand, the Si^^rn/ Sg^rn flux ratio does not 
depend on the mass (luminosity) of galaxies, analogously 
to the specific SFR or SFR per unit mass. Therefore, a 
census of star-formation activity can be performed using 
the Sis^rn/Sg^m flux ratio. 

In Fig. 5, we show normalized histograms of 
Sib^m/ Sg^rn flux ratios. The red solid line is for 
cluster member galaxies, which were selected in Section 
4.2 as galaxies in 5{z' — NA) < 0.2 from the best-fit color 
of the red-sequence, and within 1.5 Mpc from the cluster 
center. All histograms are normalized so that fractions 
add up to unity. The blue dashed line is for field galaxies 
in the same redshift range (0.9 < z < 1.7) selected using 
the photometric redshift computed by Takagi et al. (2008) 
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Fig. 1. Cluster candidates at 0.9 <z < 1.7 found in the NEP deep field. The left panels are Subaru z^-band images. The middle 
panels are AKARI A^4-band images. The scale is indicated in the left panels. Galaxies marked with red circles (2" radius) have 
consistent z' — NA color to be on the red-sequence at the detection redshift. The right panels are N4 vs z' — N4 color-magnitude 
diagram, where large black dots are cluster member galaxies. The orange dots show model galaxy colors with evolution as a function 
of redshift. The red solid line shows the color-magnitude relation of the Coma cluster at the detection redshift. The green triangles 
are galaxies which have photometric redshifts consistent to be at a cluster redshift. For reference, the small gray dots in the 
background are all galaxies within 1.5 Mpc from the cluster centre. 
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Fig. 1. (Continued). 




Fig. 1. (Continued). 
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Fig. 1. (Continued). 
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Fig. 1. (Continued). 



and Wada et al. (2008). Briefly, we used the hyperz code 
with GALEX FUV,NUV,B,V,RJ,z,J, and K as inputs 
whenever available. The resulting photometric redshift 
estimates agree reasonably well with a dozen of galaxies 
with spectroscopic redshift in the NEP fleld (Takagi et al. 
2008). We also show galaxies with photometric redshift 
consistent to be at a cluster redshift in Fig.l with the 
green triangles. Some photo-z selected galaxies are redder 
than the red-sequence. These may be dusty-galaxies or 
errors in photometric redshift. Identities of these need to 
be verifled spectroscopically. The purple dotted-dash line 
is photo-z selected cluster galaxies which have consistent 
photometric redshits with that of the cluster itself and 
within 1.5 Mpc of radius. 

Compared with field galaxies at 0.9 < z < 1.7, the 
Sih^m/ S^^rn distribution of cluster red-sequence galaxies 
are skewed toward a lower Si^^rn/ Sg^rn ratio. The median 
Si^^m/ Sg^rn ratios are 2.3 and 1.9 for the field and cluster 
red-sequence galaxies, respectively. The photo-z selected 
cluster members also have lower ratio of Si^^rn/ Sg^rn than 
the field galaxies. The median Si^^rn/ Sg^rn ratio is also 
lower with 1.4 compared with 2.3 for the field galaxies. 
According to a Kolomogorov-Smirnov test, the photo-z 
and red-sequence cluster members are different from the 
photo-z selected field galaxies with 73 and 88% of sig- 
nificance. The difference indicates that both cluster red- 
sequence and photo-z galaxies have lower star formation 
activity even at 0.9 < 2; < 1.7, quite similarly to red, passive 



galaxy clusters in the local Universe. At the same time, 
the difference in the Si^^rn/ Sg^rn ratio suggests that our 
clusters are not a coincidence of an alignment of galaxies, 
but a real clustering of galaxies in redshift space. 

5. Discussion 

Caution must be taken when the number density of our 
clusters is statistically discussed. Since we sample galax- 
ies to the detection limit of the AKARI NA magnitude, 
inevitably a narrower luminosity range (by about half a 
magnitude) is sampled for higher redshift clusters, i.e., our 
sample of clusters are not a luminosity nor mass selected 
sample although our cluster detection is machine-based, 
and can be corrected for selection effects in principle by 
using a numerical simulation. Rather we have attempted 
to the best of our data to detect as low mass clusters as 
possible toward higher redshift where optical searches are 
difficult. 

The total volume surveyed in the 0.4 deg^ NEP field at 

0. 9 < z < 1.7 is 9.1xl0^h-^Mpc^. In the local universe, 
the number density of massive clusters (~ lO^^M©) is 
3xlO"^h^Mpc-^ (Rines et al. 2007; Rines et al. 2008). At 
2: > 1, the cluster number density is about an order smaller 
than this, assuming the Press- Schechter formulation and 
the ACDM cosmology (Bahcah & Fan 1998; Carlstrom et 
al. 2002), becoming close to 3x 10"^h^Mpc-^. In Table 

1, we show 16 promising cluster candidates. We have 
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Fig. 3. The left panel shows the z^-band image of a known galaxy cluster RXJ175719. 4+663131. The scale is indicated in the 
figure. The right panel shows the A^4-band image of the same region. Galaxies marked with red circle have consistent z' — color 
to be at z=0.6909. 



searched for a volume that could include a few massive 
clusters (- 10^^ Mq) in the AKARI NEP deep field. For 
a comparison, van Breukelen et al. (2006) also found 13 
clusters in the 0.5 deg^ region of the UKIDSS survey. 

We have cross-correlated our cluster candidates with 
the positions of ROSAT X-ray detections (Henry et al. 
2006). However, there was no match except for the 
RXJ175719.4+6633 at z=0.6909. The match proves the 
validity of our cluster detection algorithm. On the other 
hand, none of our z > 0.9 clusters were matched with the 
known X-ray sources, as is expected since the flux-limit of 
the ROSAT is - 2 x lO"^"^ ergs cm-^s"^ in the 0.5-2.0 keV 
band, which requres large X-ray luminosity of > 2 x 10^^ 
erg s~^ at z > 1 to be detected with the ROSAT. The ex- 
pected X-ray luminosity of our clusters are much smaller 
than this at around 3 x 10^^ erg s~^ assuming the mass of 
< 10^"^ Mo (Maughan 2007). 

In the literature, there have been a few attempts to 
locate galaxy clusters at z > 1 using the optical-near in- 
frared data. McCarthy et al. (2007) found a compact 
group of massive red galaxies at ^ = 1.51 in 120 arcmin^ of 



the Gemini Deep Deep survey. The UKIDSS Ultra Deep 
Survey detected 13 cluster candidates with 0.61 <z< 1.39 
in a 0.5 deg^ survey (Yamada et al. 2005; van Breukelen et 
al. 2006). Castellano et al. (2007) found an overdensity of 
galaxies at z=1.6, presumably a poor cluster in the process 
of formation, in the GOODS-MUSIC catalog using pho- 
tometric redshifts only. However, these surveys still suffer 
from relatively smaller fleld of view of the near-infrared de- 
tectors of ground-based telescopes. A major step forward 
will perhaps be brought by infrared space telescopes such 
as the Spitzer space telescope and our AKARI satellite, 
which can survey large area faster without suffering from 
the atmospheric background in the infrared. Eisenhardt 
et al. (2008) have identified 106 galaxy cluster candidates 
at z > 1 using a 4.5 fim selected sample of objects from 
a 7.25 deg^ region in the Spitzer Infrared Array Camera 
(IRAC) Shallow Survey and the NOAO Deep Wide Survey 
area. They have spectroscopic confirmation for 12 clusters 
out of the 106. Although our AKARI NEP DEEP survey 
is 0.4 deg^, cluster finding in the AKARI NEP wide sur- 
vey in 6 deg^ (Matsuhara et al. 2006; Matsuhara et al. 
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Fig. 5. Normalized distribution of S15/S9 flux ratio. All his- 
tograms are normalized so that fractions add up to unity. The 
red solid line is for red-sequence cluster galaxies in the red- 
shift range of 0.9 < z < 1.7. The blue dashed line is for fleld 
galaxies in 0.9 <z< 1.7 selected using the photo-z. The purple 
dash-dotted line is for cluster galaxies in 0.9 < 2; < 1.7 selected 
using the photo-z. The median Si^/j^m/ Sg^rn ratios are 1.9, 
1.4 and 2.3 for the cluster red-sequence, photo-z cluster, and 
fleld galaxies, respectively. 

2007; Lee et al. 2007) is underway with somewhat shal- 
lower detection limit. In addition, it is our advantage that 
the AKARI fields have photometry in 9 consecutive pass- 
bands in 2-24/im {N2, N3, N4, S7, S9W, Sll, L18W 
and L24), allowing us to immediately investigate near-, 
mid-infrared properties of cluster galaxies. 

Our detection of no enhanced PAH emission in terms of 
Sib^m/ S^^m emission in 0.9 < z < 1.7 clusters is quite in- 
teresting. Cooper et al. (2006) found that among blue 
galaxies, more luminous galaxies are in the dense re- 
gions at z ^ 1. Elbaz et al. (2007) showed that the 
star formation-density relation was reversed at z^l: the 
average SFR of an individual galaxy increased with lo- 
cal galaxy density at z~l. Recently, Koyama et al. 
(2008) found possible dusty-starforming galaxies detected 
in 15/im in a galaxy cluster at z ~ 0.8. Serjeant et al. 
(2008) showed that at least a third of the Spitzer-identified 
submm galaxies at 1 < 2: < 1.5 appear to reside in over- 
densities, supporting the high-redshift reversal of the local 
starformation-density relation. 

These results indicate that galaxies in the dense re- 
gions are rapidly evolving at z~l. In the local Universe, 
the galaxy SFR is lower at dense regions (Gomez et al. 
2003; Tanaka et al. 2004). According to the downsiz- 
ing galaxy formation scenario, where stellar growth sys- 
tematically moves to lower mass systems (Cowie et al. 
1996), more massive galaxies are expected to be star 
formating at higher redshift(Kodama et al. 2004; Treu 
et al. 2005; Cimatti et al. 2006; Neistein et al. 2006). 
In this case, one might expect galaxies in dense regions 



have stronger PAH emission at z ^1. Our finding, how- 
ever, is opposite; cluster red-sequence and the photo-z 
selected cluster galaxies are already more passive even 
at 0.9 < z < 1.7. The formation epoch of these passive 
red-sequence galaxies must be at least a few Gyrs fur- 
ther in the past from 0.9 < z < 1.7. At the same time, 
our results suggest that galaxy formation epoch is quite 
earlier in the dense regions than in the field region, i.e., 
our results require even stronger downsizing galaxy for- 
mation (Cowie et al. 1996), and possibly implying the 
presence of separate physical mechanisms responsible in 
the dense/sparsely populated regions (Goto et al. 2003b). 

6. Summary 

The deep 7V4 image with the AKARI, the first Japanese 
satellite dedicated to infrared astronomy, combined with 
the 2;'-band imaging with the Subaru telescope provides 
us with an opportunity to find massive galaxy clusters at 
0.9 < 2; < 1.7, where previous searches for galaxy clusters 
based on optical photometry were difficult. By applying 
an efficient galaxy cluster finder we developed, we have 
selected 16 promising galaxy cluster candidates at 0.9 < 
z < 1.7 in the AKARI NEP deep field. 

Our cluster candidates show an obvious over-density of 
galaxies in the Subaru 2;'-band image, in addition to a 
tight color- magnitude relation in z^ — N4 vs. 7V4 (Fig.l). 
And thus, they are very likely to be a true clustering of 
galaxies at 0.9 < z <1.7. 

Using these high redshift galaxy clusters, we have in- 
vestigated the distribution of the Si^^rn/ S^^rn flux ratio, 
which is an extinction-free indicator of star-formation ac- 
tivity at 0.9 < z < 1.7. Compared with photo-z selected 
fleld galaxies in the same redshift range, both the red- 
sequence and photo-z selected cluster member galaxies 
statistically show a lower Si^^rn/ S^^rn flux ratio. This re- 
sult indicates that cluster member galaxies have less star- 
formation activity even at 0.9 < z < 1.7, suggesting that 
the formation epoch of these cluster galaxies are at higher 
redshift. 
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